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ABSTRACT 

Cellular oxidative and electrophilic stress triggers a 
protective response in mammals regulated by 
NRF2 (nuclear factor (erythroid-derived) 2-like; 
NFE2L2) binding to deoxyribonucleic acid- 
regulatory sequences near stress-responsive 
genes. Studies using Nrf2-deficient mice suggest 
that hundreds of genes may be regulated by 
NRF2. To identify human NRF2-regulated genes, 
we conducted chromatin immunoprecipitation 
(ChlP)-sequencing experiments in lymphoid cells 
treated with the dietary isothiocyanate, 
sulforaphane (SFN) and carried out follow-up biolo- 
gical experiments on candidates. We found 242 
high confidence, NRF2-bound genomic regions 
and 96% of these regions contained NRF2- 
regulatory sequence motifs. The majority of 
binding sites were near potential novel members 
of the NRF2 pathway. Validation of selected candi- 
date genes using parallel ChIP techniques and in 
NRF2-silenced cell lines indicated that the expres- 
sion of about two-thirds of the candidates are likely 
to be directly NRF2-dependent including retinoid X 
receptor alpha {RXRA). NRF2 regulation of RXRA 
has implications for response to retinoid treat- 
ments and adipogenesis. In mouse, 3T3-L1 cells' 
SFN treatment affected Rxra expression early in 
adipogenesis, and knockdown of Nrf2-delayed 
Rxra expression, both leading to impaired 
adipogenesis. 



INTRODUCTION 

NRF2 [nuclear factor (erythroid-derived 2)-like 2; 
NFE2L2] is a key transcriptional mediator of antioxi- 
dant and Phase II detoxification genes, which collect- 
ively protect the cellular milieu from oxidative and 
electrophilic stress (1,2). In response to stimuli, NRF2 
dissociates from KEAP1, translocates to the nucleus 
and dimerizes with other bZIP proteins such as small 
Maf proteins to form a transactivation complex that 
binds to antioxidant response elements (AREs) (3). 
/V>/2-deficient mice are highly susceptible to 
chemical-induced toxicity, carcinogenesis and oxidative 
burden (4), particularly in the lung, and this likely to be 
due to the loss of expression of detoxification, xeno- 
biotic metabolism and redox genes (1,2,5-9). Of 
growing importance, other gene classes including those 
involved in protein transport, ubiquitination, phosphor- 
ylation, cell cycle, growth and apoptosis have been 
identified as potentially /V>/2-dependent and are 
altered in studies with phenolic antioxidant and isothio- 
cyanate treatments (10-14). 

To identify genes that are directly regulated by NRF2 
requires evidence of binding at AREs and transactivation 
of transcription. Recent studies have used bioinformatics- 
based motif recognition (3,15-17), sequence conservation 
and gene expression information to identify genome-wide 
targets for NRF2 (18,19). Presumption of NRF2 occu- 
pancy of AREs has frequently been determined using elec- 
trophoretic mobility shift assays (EMSAs), luciferase 
constructs and immobilized oligonucleotide-binding 
assays (20), but chromatin immunoprecipitation (ChIP) 
is more definitive. Recently, Malhotra et al. (21) used 
deep sequencing of ChIP deoxyribonucleic acid (DNA) 
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(ChlP-seq) to identify several hundred new Nrf2 targets in 
mouse embryonic fibroblasts using a genetic knockout 
strategy. Although there are more than 30 human genes 
demonstrated to be NRF2 dependent, less than half of the 
binding sites (AREs) have been validated by NRF2 ChlP. 
Thus, our goals for this work are to: (i) use ChlP-seq in 
human cells to localize and validate NRF2 binding at 
both known and 'putative' binding sites (i.e. targets 
lacking conclusive binding evidence in the human 
genome); (ii) discover novel candidate NRF2-mediated 
genes; and (hi) link these new candidates to biological 
processes in new ways. 

Lymphoblastoid cell lines (LCLs) have been used to 
model exposure in human populations because of the 
rich genetic information available (22-24), and we have 
used them to examine NRF2 activation after treatment 
with the dietary isothiocyanate, sulforaphane (SFN). 
Dietary antioxidant supplementation with SFN in 
humans induces antioxidant genes in vivo (25), alters 
HDAC activity in peripheral blood mononuclear cells 
(26) and protects cells from diverse types of toxic envir- 
onmental exposures (27). For example, dietary antioxi- 
dants seem to increase resistance to oxidative stress 
and mutagenicity in human lymphocytes (28,29), and 
this is likely due to activation of the NRF2 antioxidant 
pathway. Thus, activation of the NRF2-mediated 
pathway is increasingly proposed as a way to prevent or 
treat disease, including cancer (30), cardiovascular (31), 
obesity (32,33) and neurodegenerative disease (34). In 
addition, a Phase II therapeutic trial of an NRF2 
activator has proven to be effective in chronic kidney 
disease (35). 

To identify novel genes regulated by NRF2, we 
activated NRF2 to induce binding to DNA and 
compared the results of genome-wide techniques for 
detection of NRF2 occupancy [ChlP-on-DNA micro- 
array (ChlP-on-chip) and ChlP-seq] with expression 
microarray data. This information was integrated with 
bioinformatics and independent functional analyses 
[ChlP-polymerase chain reaction (PCR) and short 
interfering RNA (siRNA) knockdown] of candidate 
gene expression in multiple cell types to support the se- 
lection of biologically plausible NRF2-dependent candi- 
date genes for follow-up. In addition to validating 
numerous previously known and putative ARE sites in 
human cells, we have characterized novel NRF2-binding 
sites near genes that are involved in processes other than 
oxidant neutralization and xenobiotic metabolism, such 
as heme metabolism, apoptosis and cell death, immune 
response and retinoid signaling by the nuclear retinoid 
X receptor alpha (RXRA). Two identified NRF2 
targets, RXRA and PPARGC1B, play central roles in 
lipid metabolism and adipocyte differentiation, and we 
further examined these genes in the context of the 
mouse 3T3-L1 adipocyte differentiation system. SFN 
treatment and Nrf2 knockdown alter expression of 
Rxra in differentiating adipocytes, and this was 
associated with impaired adipogenesis. 



MATERIALS AND METHODS 

Cell culture 

We grew human lymphoblastoid cells (Coriell) as recom- 
mended by Coriell. BEAS-2B, A549 and 3T3-L1 cells were 
grown as recommended by ATCC. Details of cell culture 
experiments and sample preparation are provided in the 
Supplementary Materials and Methods section. 3T3-L1 
cells were grown until confluent (designated as Day 0), and 
differentiation process was initiated by replacing growth 
medium with differentiation medium containing 
dexamethasone, 3-isobutyl-l-methylxanthine and insulin 
(DMI) cocktail as per published methods and stained with 
Oil Red O (32,36,37). 3T3-L1 cells with stable shRNA 
knockdown for Nrf2, Keapl and scramble shRNA were 
generated as described in Pi et al. (32). Insulin solution 
(human, 19278), 3-isobutyl-l-methylxanthine (IBMX, 
17018) and dexamethasone (D1756) were obtained from 
Sigma (St Louis, MO, USA) and d,l-SFN from 
Calbiochem. 



ChlP and ChlP-Seq 

For each tested LCL (total 6), we performed ChlP in 
biological duplicate or triplicate and combined replicates 
after validating SFN-induced NRF2 binding at known 
ARE loci for every experiment (Figure IB and Supple- 
mentary Figure S2). Details are in the Supplementary 
Methods. ChlP DNA was quantified using Picogreen 
(Invitrogen), and sequencing libraries were prepared 
using the standard Illumina ChlP-seq protocol. The 
National Center for Genome Resources (Santa Fe, NM, 
USA) created the libraries and sequenced the immunopre- 
cipitated samples on the Illumina Genome Analyzer II. 
Reads were mapped to NCBI human reference genome 
(build 36.3, March 2008) using Burrows-Wheeler 
Alignment (BWA) Tool (38). Sequencing data has been 
deposited to the National Center for Biotechnology 
Information GEO database (http://www.ncbi.nlm.nih. 
gov/geo/) under accession number GSE37589. 



ChlP-seq data analysis 

The uniquely mapped short reads were used to identify 
regions of the genome with significant enrichment in 
NRF2-associated DNA sequences, hereafter referred 
to as 'peak regions' because of their appearance 
in genome-wide density plots. The peak detection was 
performed by QuEST 2.4 software (39) using the 
'Transcription factor binding site' setting (bandwidth of 
30 bp, region size of 300 bp) and the 'stringent peak 
calling' parameters (corresponding to 50-fold ChlP to 
input enrichment for seeding the regions and 3-fold 
ChlP enrichment for extending the regions). The three 
sequenced lanes were examined individually for reprodu- 
cibility, and then unique reads were combined and 
reanalyzed to produce composite binding regions. 
ENCODE data (40) were downloaded from the UCSC 
Browser. 
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Figure 1. SFN activates NRF2 nuclear translocation, downstream ARE 
gene expression and binding to ARE genes in lymphoblastoid cell lines. 
(A) NRF2 protein levels increase in nucleus after 4-h 10-u.M SFN treat- 
ment, and these levels persist through 24 h. TBP, TATA-binding protein, 
HSP90. (B) Gene expression of known NRF2-regulated genes, HMOX1 
and NQOl, significantly increase after 4-, 8- and 24-h 10 -uM SFN 
exposures. Gene expression values are normalized to GAPDH. Bars rep- 
resent average of six independent experiments ± SEM; *P < 0.05. 
(C) ChlP-PCR analysis demonstrates that 5-h 10 -uM SFN treatment 
increases NRF2 binding at known ARE loci in the promoter regions of 
HMOX1 and (D) NQOl. A small, but significant, amount of binding was 
also seen in vehicle control-treated cells. Bars represent amplified target 
loci normalized to input DNA, as described in 'Materials and Methods' 
section. Controls: immunoprecipitation (IP) with non-specific, 
species-matched IgG antibody and amplification of loci lacking a func- 
tional ARE motif (NQOl promoter 2-kb upstream of functional AREs 
and p21 promoter region). Data in (C) are means ± SEM of seven inde- 
pendent experiments in a single cell line; (D) displays mean ± SEM of six 
independent ChIP experiments in six different LCLs. Enrichment graphs 
for several individual experiments are shown in Supplementary File 1, 
Supplementary Figure S2. Each ChIP experiment was tested in this way 
before sequencing. 



Gene expression micro array 

To measure whole-genome expression in SFN-treated 
LCLs, we used Illumina HumanRef-8 v. 3.0. Expression 
BeadChips. 

De novo motif discovery and identification of putative 
NRF2-binding sites 

To determine whether our ChlP-seq experiment allows 
unambiguous recovery of the DNA motif responsible for 
sequence-specific NRF2 protein binding, we applied de 
novo motif discovery to the ChlP-seq peak regions. 
Using the Gibbs motif sampler (41,42), we searched all 
regions for enriched sequence patterns without any as- 
sumption. To identify putative NRF2-binding sites, we 
used a position weight matrix method based on a 
curated list of known AREs, as described previously (19). 

Real-time quantitative PCR 

For ChlP-PCR, we isolated ChIP DNA, as described 
earlier and tested for enrichment of ARE loci in 
HMOX1, NQOl and other loci using the Power SYBR 
PCR Master Mix (ABI). For gene-specific expression 
analysis, we isolated total RNA isolated using the 
RNeasy kit (Qiagen) and treated with DNase. We 
reverse transcribed complementary DNA using 
Superscript First-Strand Synthesis System (Invitrogen) 
with poly-T primers following manufacturer's instructions 
and used TaqMan assays or exon spanning primers 
(Supplementary Table S6). 

NRF2 gene silencing 

We silenced NRF2 in BEAS-2B and A549 cells using a 
reverse transfection protocol using Ambion siRNAs. Nrf2 
and Keapl were silenced in 3T3-L1 cells using short 
hairpin RNA (shRNA) as described in Pi et al. (32). 
Additional details are available in Supplementary 
Materials and Methods section. 



RESULTS 

Activating NRF2 in LCLs 

After exposure to oxidative or electrophilic stress, NRF2 
dissociates from the cytosolic inhibitor KEAPl. It subse- 
quently translocates to the nucleus (Figure 1A), forming a 
heterodimer with small MAF proteins (i.e. MAFG) and 
targeting ARE-binding sequences to activate transcrip- 
tion. To identify NRF2 occupancy at novel sites, we 
used a non-cytotoxic dose of the dietary isothiocyanate 
SFN to robustly activate NRF2, induce DNA binding 
and the transcription of downstream target genes in 
human LCLs from the HapMap CEU collection 
(Figure 1A-C; Supplementary Figure SI) (6,9). 
Figure 1A displays NRF2 protein translocation to the 
nuclear cellular fraction after 10 uM SFN treatment at 4, 
8 and 24 h, and these levels were sustained for at least 24 h. 
Using real-time PCR, we verified NRF2-mediated gene 
transcription of well-known ARE genes, HMOX1 (43) 
and NQOl (15) at 4, 8 and 24 h (Figure IB). 
Importantly, using ChIP real-time PCR, we verified 
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NRF2 binding at known NRF2-binding loci in the pro- 
moters of either HMOX1 (Figure 1C) or NQOl 
(Figure ID and Supplementary Figure S2) after SFN 
treatment in each of the independent ChIP experiments 
carried out in replicate in multiple LCLs. We observed 
that in untreated cells, low levels of NRF2 protein were 
found in the nucleus consistent with NRF2 occupancy at 
the HMOX1 and NQOl promoters (Figure 1A, C and D). 
Both of these phenomena may be due to background 
oxidative stress in cells cultured at ambient oxygen level 
(44). These experiments established NRF2 activation, 
promoter occupancy and inducible downstream gene 
expression from 4 to 24 h SFN treatment in our cell 
culture model system. 

Genomic distribution of NRF2 binding 

We selected six unrelated HapMap CEU LCLs to treat 
with SFN, and following verification of NRF2 occupancy 
by ChlP-PCR (Figure ID), genome-wide NRF2 binding 
was determined with sequencing. We isolated immunopre- 
cipitated chromatin using a monoclonal antibody specific 
for human NRF2 from replicate cultures for each cell line 
treated with vehicle only (0.1% dimethylsulphoxide (also 
referred to as No Treatment, NT)) or 10 uM SFN for 5 h. 
We sequenced ChIP DNA with an Illumina Genome 
Analyzer II and aligned the 36-base sequence reads to 
NCBI human genome build 36.3 using the BWA 
program (Supplementary Table SI) (45). In addition, we 
sequenced input DNA and DNA immunoprecipitated 
with immunoglobulin G (IgG) antibody. To detect 
ChlP-seq peaks, we used the QuEST program (39) with 
its 'stringent peak calling 1 parameters and compared peak 
calling using a variety of reference sequences from input 
and IgG ChIP (Supplementary Table S2). Supplementary 
Figure S3 displays an example of how uniquely mapped 
sequence reads upstream of KEAP1 were distributed on 
both forward and reverse strands and then converted to 
density profiles for peak calling. We examined significant 
peaks in individual sequencing runs and only considered 
further those peak regions appearing in at least two 
sequencing runs from independent ChIP experiments 
(Supplementary Figure S4 and Supplementary Data set 
1). To increase the sequencing depth and to further 
assess peak quality, we pooled the uniquely mapped 
reads from the three SFN and NT sequencing runs 
(Supplementary Figure S4). We identified 849 NRF2- 
occupied regions after SFN treatment (mean reads/ 
peak = 22.7; total reads within peaks = 19 345) and for 
untreated, vehicle-only samples, 205 peak regions (mean 
reads/peak = 8.9; total reads within peaks = 1843). Of the 
205 peak regions in the vehicle-only treated samples, 39 
(19%) overlapped with peak regions in the SFN-treated 
samples. The 4.1 -fold increase in the number of peak 
regions and 10.4-fold increase in uniquely mapped reads 
within peaks for the SFN-treated versus untreated cells 
were consistent with both the increase of NRF2 in the 
nucleus and the increased expression of known 
NRF2-mediated genes in SFN-treated cells. Additionally, 
many peak regions in both samples were just above back- 
ground levels, therefore we also assessed the number and 



frequency distribution of ChlP-seq reads per peak in the 
SFN-treated sample and considered peaks with more than 
15 unique reads (upper 95% of the mean) as high confi- 
dence (HC) peaks (242 of 849 peak regions). Only 19 of 
205 peak regions in the untreated sample had more than 
15 unique reads. 

We assessed the location of sequence reads and ChlP-seq 
regions in relationship to the nearest gene transcriptional 
start site (TSS) based on the NCBI human RefSeq gene 
database (26 733 genes), Supplementary Figure S5A-S5C. 
NRF2-binding sites have been reported as ri.v-acting 
elements located at an average distance of ~1800 bp from 
the gene TSS, based on 39 currently known functional 
human AREs (Supplementary Table S3). Among the 
SFN-treated ChlP-seq peaks, we observed that 716 of the 
849 (84%) peak regions were within 100 kb of a gene TSS 
(Supplementary Figure S5B). Of these regions, 37% (268) 
of these regions were within 5 kb of gene TSS, and this 
represents a ~ 10-fold enrichment of peaks near TSS over 
background (37% versus 3.6% of the genome maps to 
within 5kb of a TSS). Similarly, when we consider the 
HC peak regions, we observed 205 of 242 (85%) were 
within 100 kb of TSS, and 28% were within 5kb of the 
TSS. When we mapped NRF2 ChlP-seq peaks relative to 
gene structural features, most peaks (64%) were located 
within 10-kb upstream to 10-kb downstream of genes 
(Supplementary Figure S5C). Overall, these data suggest 
that most NRF2-binding sites, including the known 
AREs, are indeed cis acting, but there are certainly many 
HC-binding sites (37/242, 15%) that are quite distant from 
aTSS(>100kb). 

Binding motif analysis and overlap with NF-E2 binding 

To determine whether our ChlP-seq experiment recovered 
the DNA motif assumed to be responsible for sequence- 
specific NRF2 protein-DNA binding (Figure 2A), we 
applied de novo motif discovery to the ChlP-seq peak 
regions. Using the Gibbs motif sampler (41,42), we 
searched all regions for enriched sequence patterns without 
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Figure 2. ChlP-seq peaks contain NRF2-binding motif. (A) NRF2- 
binding motif based on known AREs and (B) most enriched 
sequence pattern in the 849 ChlP-seq peak regions based on a 
de novo analysis using Gibbs motif sampler. 



7420 Nucleic Acids Research, 2012, Vol. 40, No. 15 



any assumption. The top enriched de novo motif (Figure 2B) 
closely matches the core NRF2 consensus binding motif 
(Figure 2A), described previously in many studies 
(3,16,19). Using a position weight matrix based on known 
AREs (19), we determined that 643 of 849 (76%) peak 
regions obtained from the SFN-treated samples had at 
least one ARE. Of the 242 HC NRF2-binding regions, 
we observed that 231 peak regions (95%) contained one or 
more AREs. The lower ARE frequency (76%) in the lower 
confidence peak regions (composed of 1 5 reads or less) may 
be due to a higher number of false-positive signals among the 
lower confidence peaks. The locations of >90% of AREs 
identified were within 100 bp of the peak maximum 
(Supplementary Figure S5D). 

NRF2, encoded by the NFE2L2 gene, is closely related 
to the erythroid-specific regulator NF-E2, which has a 
similar API-like binding motif (TGA(G/C)TCA). Given 
the similarity of their binding motifs, we reasoned that 
these two DNA-binding proteins should share ChlP-seq- 
binding regions. We examined the overlap between 
NRF2-binding regions in lymphoid cells and NF-E2- 
binding regions recently determined in the erythroid cell 
line K562 by the Snyder lab in the ENCODE project (40) 
and found that 112 of 242 HC peaks (46%) were shared 
between these related proteins (5.6-fold enrichment over 
the mean overlap [8%, 20/242] between NRF2-binding 
sites and the 14 ENCODE TFs determined by ChlP-seq; 
Fisher's exact test, P< 0.0001), further supporting the 
functionality of the identified NRF2-binding regions. 
Although NF-E2 is exclusive to erythroid cells, NRF2 is 
widely expressed among tissues. Notably, NRF2 has 
recently been identified as a key regulator of oxidative 
stress response in erythroid cells in sickle cell disease (46). 

Localization of binding sites in known and novel 
NRF2-dependent genes 

Figure 3A SFN track displays an example of 
SFN-induced NRF2 ChlP-seq peak region localized 
over a known functional ARE in the NQOl promoter 
(47). Supplementary Data set 1 lists genes with significant 
ChlP-seq peaks in the SFN and NT sample. Seventeen of 
the top 25 genes ranked as the most significant binding 
peaks have been previously described as bone fide or 
putative NRF2-dependent genes, and these ChlP-seq 
data provide further evidence of NRF2-binding-mediated 
regulation for these genes (e.g. FTL, FTH1, HMOX1, 
NQOl, TXN. GCLC, TXNRD1, ME1 and PRDX1). In 
addition, we identified binding locations in several 
putative ARE genes previously reported as part of an 
NRF2-regulated gene expression signature in airway 
cells (PIR, TKT, ABCB6, SQSTM1, GSR and 
TALDOl) (48). Among human genes considered to have 
functional ARE sequences, we observed that some 
(e.g. PRDX1 and PIR) (Figure 3B and C) displayed 
NRF2-binding peaks in locations different from those 
described in the literature. PRDX1 displayed a prominent 
peak over a consensus matching, evolutionarily conserved 
ARE sequence located in the proximal promoter (PRDX1, 
chrl:45 760275) and PIR displayed a strong peak in the 
proximal promoter (Figure 3C, chrX: 15 421 342). In both 
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Figure 3. Examples of NRF2 ChlP-seq peaks at AREs in the gene 
promoter regions of (A) NQOl, (B) PIR and (C) PRDX1. Control 
reads from IgG IP are shown in the first track (IgG), vehicle-treated 
samples in the second track (Veh) and SFN-treated samples in the third 
track (SFN). Ticks represent reported ARE sites, based on previous 
literature (Supplementary File 1 and Supplementary Table S4) or 
observed AREs. 



cases, we did not detect binding at the locations reported 
by other investigators (48,49). Importantly, among the 
most significant HC peak regions (Supplementary Data 
set 1, SFN), there are many regions that have never 
been functionally linked with NRF2. These novel peak 
regions are nearby the TSS of genes, which include 
RXRA, AKIRIN2, GPI, KIFC3, ANKRD11, SH3TC1, 
CPEB2, CPEB3, PPARGC1B, UNKL and BMP10. 
Interestingly, the location of the peak region in the 
RXRA (RXRA) gene suggests that an alternative TSS 
may be used in LCLs. RXRA has several predicted alter- 
native start sites and several reported alternatively spliced 
messenger ribonucleic acids (mRNAs). We observed a 
large NRF2-binding peak located over a conserved ARE 
sequence (Figure 4A-C) located in intron 2 near the TSS 
of one of the shorter transcripts. In Figure 4A and in the 
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lower tracks on Supplementary Figure S6A-S6D, we 
display ENCODE ChlP-seq data for histone H3K4me3 
(active promoter) and RNAPol2 (transcription) for 
GM12878 (LCL) and HepG2 (hepatocarcinoma). The 
relative locations of these signals, near 5' TSS for 
HepG2 and near alternate TSS in GM 12878, suggest 
that transcription of RXRA may start at different TSS 
in these two lines. 

We did additional ChIP experiments to prioritize 
candidates for further analysis, including a ChlP-on-chip 
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experiment using the Agilent Human Promoter array, and 
determined the overlap with ChlP-seq peaks (Supple- 
mentary Data set 1). Three examples of candidate NRF2 
genes (OSGIN1, HTATIP2 and FECH) with highly sig- 
nificant ChlP-seq peaks, ChlP-on-chip signals and good 
ARE sequence motifs within the ChlP-seq peaks are 
shown in Figure 5A-C. For each of the binding peaks 
associated with the eight genes shown in Figures 3-5 
and Supplementary Figure S3, as well as eight additional 
genes listed in Supplementary Table S6, we further 
validated binding by ChlP-PCR in independent experi- 
ments. We detected NRF2 occupancy at 10x10 HC 
peaks and 5x5 lower ranked peaks (S6). Due to the 
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Figure 4. ChlP-seq peak region near an alternative TSS for the RXRA 
gene. (A) ChlP-seq peak region is mapped ~80-kb downstream of the 
RXRA TSS as defined by NCBI genome build 36.3. Tracks displaying 
ENCODE H3K4me3 and RNAPol2 ChlPseq data indicate open chro- 
matin and active transcription near the alternate start site 
(Supplementary Figure S6). The peak region is within 1 kb of an alter- 
native TSS of RXRA, seen in more detail in (B). (C) Sequence under 
the peak region contains a highly conserved ARE sequence, which 
closely matches the consensus sequence (blue) as shown by the 
sequence logo. Mouse and rat nucleotides not matching the human 
sequence are displayed in red. 
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Figure 5. Peak regions that matched in ChlP-seq and ChlP-on-chip. 
Three examples of genes, OSGIN1 (A), HTATIP2 (B) and FECH (C), 
which have not been previously described as regulated by NRF2. Insets 
display ARE motif under the ChlP-peak regions. Lower case nucleo- 
tides represent mismatch with ARE consensus. 
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importance of RXRA in numerous physiological 
processes, we carried out experiments to validate NRF2 
binding and inducible expression in additional cell types. 
Significant NRF2 binding in RXRA was detected in un- 
treated LCLs (Figure 4A and B; Supplementary Figure 
S6E), in SFN-exposed BEAS-2B bronchial epithelial 
cells (Supplementary Figure S6E) and mouse 3T3-L1 
cells (Supplementary Figure S6H). In addition, using 
ENCODE data for the GM 12878 cell line, we detected 
nucleosome displacement indicated by collocation of 
H3K4mel displacement and DNase hypersensitivity at 
the NRF2-binding site (Supplementary Figure S6F). 

SFN-mediated gene expression 

We used Illumina human Ref-8 microarrays to assess 
NRF2-mediated gene expression in the six sequenced 
lines and in the 60 unrelated CEU HapMap LCLs after 
24 h exposure to 10 uM SFN as part of a related, ongoing 
genetic study in our laboratory (50). Comparing mean log 2 
expression values using f-tests and fold-change values for 
treated versus untreated across all cell lines, we found 508 
genes changed by 1.3-fold or greater with false discovery 
rate (FDR) (/-value < 0.05 (Supplementary Data set 1). 
Similar results were obtained when analysis was carried 
out using only the six cell lines that were sequenced but 
with less statistical confidence (e.g. n = 6 versus n = 60). 
Among genes with both ChlP-seq peaks and gene expres- 
sion changes, there were significantly more ChlP-seq peak 
regions near up-regulated genes (20.6%; 60/291) than 
down-regulated genes (4.6%; 10/217; P< 0.0001, Fisher's 
exact test). None of the down-regulated genes displayed 
HC ChlP-seq peaks. These data support findings in the 
existing literature suggesting that NRF2 is primarily a 
positive regulator of gene transcription rather than a 
suppressor (51). Considering the 242 HC peak regions, 
expression of 97 of the nearby genes was near or below 
the array detection limit in the LCLs (log 2 expres- 
sion < 6.0), including some with highly significant NRF2 
binding (e.g. HBE1, KIFC3, BMP10 and ANKRD11). 

A number of known NRF2-regulated genes (FTL, 
PRDX1, GSTP1, PSMA3 and TKT) were highly expressed 
but showed minimal change following treatment (<30%). 
Based on the presence of strong H3K4me3 marks (enrich- 
ment at enhancer/promoter elements, ENCODE data not 
shown), expression of these genes may already be at 
maximal levels under NT conditions. They were among 
the top 2% highest expression values in the untreated 
samples, and after SFN treatment, despite the very large 
increase in NRF2 binding at their AREs (5- to 20-fold 
increase), only a small increase in expression signal was 
observed. Presumably, there are other determinants, such 
as tissue-specific cofactors, negative feedback loops, epi- 
genetic or signaling mechanisms, which affect both basal 
expression and NRF2-mediated transcriptional regulation 
of these highly expressed genes in lymphoblastoid cells. 
Among all genes with expression and ChlP-seq peaks, 67 
showed both NRF2 occupancy and SFN-altered expres- 
sion >30%, and we carried out further analysis of this 
group. RXRA induction by SFN was validated in a time 
course experiment in LCLs (Supplementary Figure S6I) 



and in several additional cell lines (HepG2, Supplementary 
Figure S6I; 3T3-L1; Figure 6D and F). 

MicroRNAs 

We identified HC ChlP-Seq peak regions in the vicinity of 
several microRNAs (miR-365-l/miR-193B, miR-617, 
miR-592, miR-29B-l, miR-1207, miR-32, miR-181C, 
miR-200C and miR-550-1), the most significant peak 
occurring near the mir365-l/mirl93B cluster (Supple- 
mentary Figure S7A-S7C). Among these, only miR-29B-l 
has been linked with oxidative stress (52), although 
miR-365-1 was observed to be induced in bronchial epithe- 
lial cells of smokers compared with nonsmokers (53). Many 
of these microRNAs have been detected in tumors and 
decreased levels of miR-365/193B (54), miR-617 (55), 
miR-200C (56) and miR-29B (57) have been linked with 
tumor aggressiveness in a variety of cancers. Using 
real-time PCR assays, we tested expression of miR-365-1/ 
miR-193B, miR-617, miR-592, miR-29B-l, miR-1207, 
mir32, miR-181C, miR-200C and miR-550-1. Of these, 
only miR-365-1, miR-193B, miR-181C and miR-29B-l 
showed appreciable expression in lymphoblast cells. 
MiR-365-1 displayed variability after exposure (trending 
up and then significantly down), but only miR-29B-l 
showed significant SFN-induced change relative to no treat- 
ment (Supplementary Figures S7D-S8C). The reduced 
expression of miR-29B- 1 in response to SFN was consistent 
with Luna et al. (52) who observed oxidant exposure 
suppressed miR-29B-l. 

Pathway analysis 

SFN is a potent NRF2 activator, but it can also regulate 
gene expression through a number of other mechanisms 
including histone modification and alteration of MAPK, 
NF-kB or AP-1 pathways (58). We did pathway analysis 
to assess whether there were gene sets (or pathways) 
enriched in the differentially expressed genes and/or genes 
with NRF2 ChlP-seq peaks. We used Ingenuity Pathway 
Analysis (IPA) to analyze gene enrichment in toxicity 
response pathways for all genes with SFN-altered gene 
expression (508 genes) and the smaller set of genes with 
both expression changes and ChlP-seq peak regions 
within lOOkb of the TSS (Supplementary Table S4). Five 
pathways were significantly enriched when examining all 
SFN-responsive genes (Supplementary Figure S9. Open 
bars. Significance threshold shown represents Fisher's 
exact test with FDR correction, P < 0.05). These were 
oxidative stress response mediated by NRF2; cell cycle: 
G2/M DNA damage checkpoint regulation; cholesterol 
biosynthesis; LXR/RXR activation and aryl hydrocarbon 
receptor signaling. These data clearly indicate that the 
NRF2 pathway is involved in SFN-mediated gene respon- 
siveness, but based on expression, several other processes 
are altered including cholesterol biosynthesis and RXR 
activation. When we limited analysis to the 67 SFN- 
responsive genes with ChlP-seq peak regions (dark bars), 
only the two pathways related to oxidative stress were 
significantly enriched in the IPA analysis. Genes with 
significant expression changes but no NRF2 binding 
at AREs may represent downstream targets of 
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Figure 6. Effects of SFN treatment on 3T3-L1 adipocytes during differentiation. (A) Table displays eight treatment conditions (SFN ± starting and 
ending on different days). Growth at day 8 is shown in photographs at two magnifications. Treatment for 2, 4, 6 or 8 days inhibits adipogenesis. 
Treatments starting at day 2 or later have less effect. (B) PparG2 expression by real-time PCR increases rapidly after day 3 under standard 
differentiation conditions indicating commitment to terminal adipocyte state. Continuous SFN treatment prevents differentiation. (C) Rxra expres- 
sion is higher under SFN treatment conditions. (D) Expression of PgclB, the PPAR coactivator (homolog of human PPARGC1B), is similar to 
Pparg2 both untreated and treated. (E) Gata3, a negative regulator of adipogenesis, dips initially, then increases dramatically by Day 3 with SFN 
treatment. (F) Rxra expression was tested in 3T3-L1 adipocytes with shRNA knockdown of Nrf2 (Nrf2-Kd, green) and Keapl (Keapl-Kd, red). 
Keapl knockdown genetically activates Nrf2. The Rxra expression profile under Keapl knockdown (red triangles) is very similar to that observed for 
SFN treatment (C, red squares). Nrf2 knockdown (green) delays the rise in Rxra expression. This is accompanied by impaired adipogenesis (32). 



NRF2-regulated genes. To examine the 67 SFN-responsive 
genes with ChlP-seq peaks in more detail, we used the 
Database for Annotation, Visualization and Integrated 
Discovery (DAVID) functional annotation clustering tool 
(59,60). Using this tool to integrate information from Gene 
Ontology, BioCarta, KEGG and more than 40 other 
annotation categories, we identified broader functional 
categories (Table 1 and Supplementary Data set 1 
DAVID). The most enriched functional categories were 
those involved in redox homeostasis and heme metabolism, 
vitamin metabolic process, oxidation reduction and 
regulation of apoptosis and cell death. 

Selection and validation of human NRF2 candidate genes 

To select promising novel NRF2-regulated gene candidates 
for follow-up under NRF2 knockdown conditions, we 
identified genes with matched peak regions in both 
ChlP-seq and ChlP-on-chip data sets that were within 
2.5 kb of the TSS and possessing an ARE within the 
peak. This subset included 29 gene candidates, 10 of 



Table 1. DAVID functional annotation clustering of 67 NRF2-bound 
genes with gene expression changes" 



Clustered category Enrichment score 



Cellular redox homeostasis 2.59 

Heme metabolism 2.39 

Vitamin metabolic process 2.37 

Oxidoreductase activity 2.34 

Apoptosis and cell death 2.33 

Iron binding (oxidoreductase activity) 1.78 

Glutathione metabolism 1.62 



"See Supplementary Data set 1 for individual categories and gene lists. 



which had known or putative functional AREs (three 
known NRF2 targets with peaks ~4.5kb from the TSS 
were also included, HMO XI, FTH1 and TFE3). An 
example of one of these candidates was FECH, the gene 
that codes for the enzyme ferrochelatase that carries out 
the final step in heme biosynthesis. FECH displayed both 
ChlP-seq and ChlP-on-chip peaks, was induced 1.3-fold 
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and had an ARE in its proximal promoter (Figure 5C). 
Given the coverage limitations of ChlP-on-chip, we also 
selected 14 genes that had a ChlP-seq peak region within 
2.5 kb of the TSS, contained an ARE motif and exhibited 
transcriptional regulation by SFN but lacked a 
ChlP-on-chip result. Among these was RXRA, which had 
a highly significant ChlP-seq peak region, contained a 
highly conserved ARE and displayed up-regulation of 
gene expression (1.5-fold). The ChlP-seq peak was not 
close to the primary NCBI-defined TSS but was within 
1 kb of an alternative TSS (Figure 4). We also chose to 
examine a pair of genes, NCAPD2 and GAPDH, with a 
binding region between them containing a highly conserved 
ARE, to test whether they 'shared' regulation by NRF2. 
These 14 additional tested genes included two well-known 
NRF2 pathway genes, MAFG and KEAP1 (KEAP1, 
Supplementary Figure S3), that had binding sites identified 
in the mouse (61,62) but not in the human genome. This set 
of genes is not intended to be comprehensive but represents 
a sampling of new candidates. 

Due to the importance of the NRF2 pathway in the 
airway response to oxidants (e.g. hyperoxia), we tested the 
impact of NRF2 silencing on expression of candidate genes 
using two airway cell lines (A549 and BEAS-2B) that are 
commonly used for NRF2-related oxidative stress studies. 
Effective silencing of NRF2 in LCLs was not achieved with 
siRNA or shRNA. However, after transient siRNA 
silencing of NRF2 transcripts, A549 showed 94% reduction 
and BEAS-2B showed 86% reduction of expression 
compared with non-specific siRNA-transfected controls. 
NRF2 silencing led to significantly reduced expression of 
six of the seven known and five of the six putative 
NRF2-regulated genes in both A549 and BEAS-2B cells 
(P < 0.05, ?-test; Supplementary Table S5). For many 
genes, we observed mRNA reduction with NRF2 silencing 
at both baseline and after SFN treatment consistent with 
high baseline NRF2 activity in these cell lines (63). NRF2 
silencing significantly altered expression of 20 of 30 of the 
candidate genes in both A549 and BEAS-2B cells including 
RXRA (Supplementary Table S5). Three genes (GPD2, 
NDUFAF4 and HIST1H4H) had reduced expression with 
NRF2 silencing in one of the lines only; four did not change 
and three were not expressed. Interestingly, we observed 
strongly increased MT2A (metallothionein 2 A) expression 
in NRF2-silenced A549 cells at baseline and a larger 
increase under SFN treatment. This suggests the possibility 
of direct negative regulation of MT2A by NRF2 in A549 
cells. This validation set is not comprehensive, and add- 
itional genes remain to be validated by silencing and 
ChlP-PCR, but these experiments reveal a number of 
potentially important new NRF2 target genes, particularly 
RXRA. Further analysis will be needed to establish their 
functional importance within the NRF2-regulated 
pathway. 

NRF2 activation and RXRA expression in adipocyte 
differentiation 

We identified strong NRF2 binding within RXRA and 
observed SFN-induced RXRA gene expression in 
multiple human cell lines following SFN treatment. In 



addition to RXRA, the PPAR coactivator PPARGC1B 
also displayed NRF2 binding, suggesting a novel role 
for NRF2 in pathways regulated by retinoids and 
involving RXRA. RXRA function is not well 
characterized in the cell types we initially investigated, 
but NRF2-binding sites in RXRA displayed 100% 
human/rodent homology (Figure 4). Thus, to explore the 
functional regulation of Rxra by Nrf2 activation, we chose 
to use the well-established mouse 3T3-L1 adipocyte differ- 
entiation model that is regulated by the Pparg/Rxra 
heterodimer. This model system has been used to impli- 
cate NRF2 in adipocyte biology (32,33,37), although the 
mechanism has remained somewhat inconclusive. To 
investigate whether Nrf2-dependent perturbation of 
Rxra could impact adipogenesis, we induced 3T3-L1 
cells to differentiate and treated the cells with SFN for 
varying times, as shown in Figure 6A. Treatment of 
3T3-L1 cells with 10 urn SFN for 2, 4, 6 and 8 days 
produced increasing levels of inhibition of differentiation 
(Figure 6A, tracks b-e). A single SFN dose on Day 1 
produced detectable inhibition (Figure 6A, track b), 
whereas repeat doses (every 2 days) produced complete 
inhibition (tracks c-e). Notably, treatments starting at 
Day 3, or after, had much less effect, indicating that the 
inhibitory effect of SFN activation of Nrf2 occurs early in 
the cascade of signals driving adipogenesis. A similar 
pattern was observed with 5 um SFN but with reduced 
effect (data not shown); toxicity was not evident. Nrf2 
ChlP-PCR of the Nrf2-binding site in the Rxra gene in 
SFN-treated 3T3-L1 cells displayed strong enrichment 
relative to NT (Supplementary Figure S6H). 

Figure 6B displays the real-time PCR profile of the 
master regulator of adipocyte differentiation, Pparg2, in 
differentiating adipocytes over 7 days under standard 
DMI conditions with and without SFN (NT vs. SFN; 
these are shown as conditions 'a' and 'e' in the Figure 6A 
table). Under standard differentiation conditions (NT), 
Pparg2 expression rapidly increases after Day 1 to a 
maximum of ~ 150-fold as the cells become committed to 
terminal differentiation. However, under SFN treatment, 
Pparg2 remains low, consistent with failing to progress to 
terminal differentiation. A similar pattern was observed for 
Cebpa (Supplementary Figure S10A). Oxidative stress 
genes known to be regulated by Nrf2 (Nqol and Gstal) 
showed modest but significant mRNA changes under 
differentiation conditions but displayed large, rapid in- 
creases during the first 3 days with SFN treatment and 
then declined to baseline levels (Supplementary Figure 
S10B-S10D). Expression of the PPARG coactivator, 
Pgclb (Figure 6E; mouse homolog of PPARGC1B) 
mirrored Pparg2, suggesting its regulation is tied to 
differentiation status rather than Nrf2 activation. 

Rxra expression, heterodimerization with Pparg, binding 
to DR1 sites and transactivation of adipocyte-specific genes 
are critical events in adipogenesis. Rxra expression was 
increased under SFN treatment relative to NT (standard 
DMI conditions) between Days 1 and 5 (Figure 6C, red 
squares and Supplementary Figure S10D). Interestingly, 
transcription factor Gata3, a negative regulator of adipo- 
cyte differentiation (64,65), was dramatically induced by 
SFN between Days 1 and 3 of treatment, and such levels 
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of Gata3 have been demonstrated to suppress adipogenesis 
through direct binding to the Pparg promoter and inter- 
action with Cebpa protein. We also used shRNA 
knockdown of Nrf2 and Keapl (the negative regulator of 
Nrf2) and examined the real-time PCR profiles of Rxra in 
differentiating adipocytes under standard DMI conditions 
(Figure 6F). Relative to Wt-Scr, Rxra expression is reduced 
in Nrf2-Kd (green) and increased in Keapl -Kd (Keapl 
knockdown produces constitutive activation of Nrf2). 
The delayed induction of Rxra expression in Nrf2-Kd 
cells may be related to the impaired adipogenesis 
observed in Nrf2-Kd cells demonstrated by Pi et al. (32). 
Importantly, the expression pattern observed for genetic 
activation of Rxra in the Keapl knockdown cells 
(Figure 6F, red triangles) was consistent with that 
observed in SFN-treated cells (Figure 6D). 



DISCUSSION 

We present the first genome-wide binding data for human 
NRF2, the transcription factor that regulates 
anti-oxidative response. Following activation of NRF2 
with the dietary antioxidant SFN, we used chromatin 
immunoprecipitation, massively parallel sequencing, 
microarray technologies, siRNA and bioinformatics to 
identify novel NRF2 human gene targets and to verify 
known or putative AREs. ChlP-seq identified 849 total 
binding regions and 242 HC regions in SFN-treated 
lymphoblastoid cells; most are novel targets in humans, 
and 110 regions overlap with ChlP-seq identified Nrf2 
target genes in mouse embryonic fibroblasts (21). Within 
these binding regions, a de novo motif analysis identified 
the most enriched motif, and it perfectly matches the core 
ARE motif. Within the HC NRF2-bound peak regions, 
96% contained one or more ARE sequences, and many, 
such as PRDX, PIR, RXRA, KEAP1, mir-29B and 
mir-365-1, were highly conserved between humans and 
rodents (Figures 3B, C and 4C, Supplementary Figures 
S3, S7C and S8B). The high frequency of AREs in these 
peaks and their position near the peak maximum support 
functionality for these binding sites. In addition, 50% of 
these peaks colocate with ChlP-seq regions determined for 
the structurally related erythroid-specific factor, NF-E2. 
The large number of genes that display NRF2 binding 
but no change in gene expression suggests that additional 
factors, perhaps coactivators and chromatin state, are 
involved in transcription regulation of this pathway. In 
addition, there are many genes that display expression 
changes but lack NRF2 binding, and these may be 
regulated indirectly by changes in the redox environment. 

RXRa, with the sixth most significant binding peak, is a 
notable addition to the NRF2-regulated pathway. It is an 
important nuclear retinoid receptor that heterodimerizes 
with RARa, RARp\ RXRp\ LXR, VDR and PPARs to 
mediate the transcriptional regulation of genes involved 
in diverse cellular processes including cell differentiation, 
cell cycle, growth and response to steroids, fatty acids and 
vitamins. Although we present the first evidence that NRF2 
regulates a retinoid receptor, our finding is not the first 
report of an interaction between the retinoid pathway and 



NRF2. Recent evidence suggests that oral retinoic acid may 
inhibit NRF2-mediated protective responses to carcinogen 
exposure (66). In addition, we and others (32,67) have 
recently identified that Pparg, the RXRA binding partner, 
is a Nrf2-regulated gene in mice. In this work, we observe 
highly significant NRF2-binding peaks near malic enzyme 
1 (ME1), a PPAR signaling pathway gene involved in 
lipogenesis and the PPAR gamma coactivator IB 
(PGC1B and PPARGC1B). In mouse studies that focus 
on dietary fat, activation of the NRF2 pathway prevents 
high-fat diet-induced obesity in mice, and although the 
mechanism is not well understood (33), a number of 
reports suggest that the effect is through the PPAR signal- 
ing pathway (32,68). In addition, NRF2 seems to be essen- 
tial for normal adipocyte differentiation in mice (32). We 
observed that activation of Nrf2 in adipocytes, either by 
SFN or genetically (Keapl-shRNA), induced Rxra 
relative to standard conditions (Figure 6B and D). In 
SFN-treated cells, after Day 1, this was accompanied by a 
rapid increase in Gata3, the negative regulator of 
adipogenesis. Considering our evidence that NRF2 activa- 
tion and binding mediate upregulation of RXRA in 
lymphocytes and adipocytes, and the clear evidence that 
PPARG/RXRA heterodimers are critical regulators of adi- 
pocyte differentiation (36,69), we hypothesize that 
NRF2-mediated regulation of RXRA expression in adipo- 
cytes may affect PPARG-regulated lipid synthesis and fat 
accumulation. If altered RXRA expression affected RXRA 
occupancy (36), this may explain how NRF2 activation 
inhibits fat accumulation (33). Recent evidence indicates 
that the ratio of RXRA to RXRB in adipocytes is import- 
ant for PPARG function (70), and perturbation of this ratio 
could be the mechanism. Thus, altered expression of RXRA 
after exposure to NRF2-activating stimuli could impact a 
number of outcomes including lipid synthesis and the 
cellular response to endogenous or therapeutic retinoids. 
A more detailed investigation of NRF2-mediated effects 
on RXRA expression during adipocyte differentiation is 
warranted. 

Our ChlP-seq analysis also identified peak regions near 
the TSS of many known or putatively identified 
NRF2-regulated genes. However, we observed the location 
of these peak regions did not always match the previously 
published positions (Supplementary File 1, Supplementary 
Table S4 and Figure 3A-C). For the PRDX1 and PIR 
ChlP-seq sites, comparative genomics analysis reveals 
consensus matching, evolutionary conserved AREs under 
the peaks (Figure 3B and C), whereas the reported NRF2- 
binding sites (48,49) have a poor match with the consensus 
and displayed little or no conservation. While not every 
human ARE displays evolutionary conservation, a prepon- 
derance of established loci do (19,71), strongly supporting 
functionality for the locations identified in this report. 

Our analysis has revealed a number of additional novel 
gene candidates that seem to expand the role of NRF2 in 
regulating oxidoreductases (AIFM2, GPD2, HTATIP2 
and NDUFAF4) and metabolism genes (GNPDA1), 
including those involved in heme metabolism (AMBP, 
FECH). It is well known that NRF2 mediates the regula- 
tion of several genes involved in iron metabolism (FTL, 
FTH1) and the heme degradation pathway (HMOX1). 
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However, NRF2-dependent regulation of the gene 
AMBP, a protein that binds heme, degrades it and is 
induced by reactive oxygen species (ROS) in erythroid 
cell lines (72), is a novel finding. In addition, we demon- 
strate for the first time that expression of FECH 
(ferrochelatase), the enzyme catalyzing the final metabolic 
step of heme biosynthesis IX, is regulated by NRF2. The 
presence of NRF2 binding in the globin locus, near HBE1 
and HBG1 is also notable, given the potential role of 
NRF2 in modulating sickle cell disease (46). 

Numerous genes identified by our ChlP-seq experiments 
suggest a role for NRF2 in cell cycle regulation and tumor 
growth. Recent demonstration that genetic activation of 
NRF2 in Keapr /_ MEFs can enhance cell proliferation 
(21) and the identification of tumor mutations that 
activate the NRF2 transcriptional pathway (73-75) are con- 
sistent with these findings. DAVID analysis of SFN- 
regulated genes that possessed ChlP-seq peak regions 
clearly indicated enrichment of functional annotation 
categories involved in cellular growth and death processes 
(Table 1). Many of the genes in this category, such as 
HTATIP2, OSGIN1 (76) and MAFG (77) have dual roles 
in oxidative stress response and cell regulation. Among 
these, HTATIP2 is an oxidoreductase that induces apoptosis 
under oxidative stress conditions by stabilization of p53 
mRNA (78). This interconnection between antioxidant 
response, redox signaling and cell cycle regulation under- 
scores the complexity and context dependence of 
NRF2-modulated responses. NRF2 activation is protective 
in countless exposure scenarios, and NRF2 activators 
inhibit growth of some tumor cell lines (14,73,74); 
however, constitutive activation in tumors by mutation 
drives proliferation and resistance to treatments. Genes, 
such as SLC3A2, implicated in teratocarcinoma formation 
in mice (79) and adenocarcinoma (80) and TFE3, regulator 
of the growth factor transforming growth factor beta, are 
further examples of proliferation and differentiation- related 
genes that are involved in tumorigenesis. Translocations in 
the TFE3 promoter region have been associated with 
sporadic cases of renal cancer (81). We also identified 
NRF2 binding near a number of microRNAs, some of 
which are implicated in tumor aggressiveness, and we 
observed NRF2-mediated suppression of miR29B following 
SFN treatment. Although follow-up of these findings is 
beyond the scope of this study, further characterization of 
the impact of NRF2 activation on genes involved in tumor 
growth is essential. 

Finally, we demonstrated NRF2 occupancy of genes 
involved in the activity and homeostasis of NRF2 itself, 
such as KEAP1, MAFG and SQSTM1. Induction of 
MAFG, the heterodimeric partner of NRF2, is mediated 
by NRF2, but binding had not been validated (17,82). The 
autophagy receptor, SQSTM1 /p62, is both a target for 
NRF2 and a positive regulator of NRF2 (83). Of particular 
interest is the observation that NRF2 binds an ARE motif 
in the KEAP1 promoter. Sustained NRF2 activity levels 
likely lead to NRF2-dependent production of KEAP1, 
which in turn leads to suppressed NRF2-mediated gene 
expression through a negative feedback mechanism. 

Although NRF2 is the master regulator of oxidative 
stress response, including during human erythroblast 



differentiation (46), it is increasingly evident that NRF2 
mediates the expression of genes involved in other diverse 
processes (21,68,84). For example, Nrf2 has been 
identified as a positive regulator of Notch 1 (85), a 
negative regulator of osteoblast differentiation in mice 
(86) and necessary for adipocyte differentiation (32). Our 
experiments have further extended the role of NRF2 sig- 
naling in adipogenesis. NRF2 regulation of RXRA in 
adipogenesis suggests the potential for NRF2 activators 
to therapeutically manipulate RXRA in numerous 
retinoid-mediated pathways. Further characterization of 
how these new targets fit into an expanding NRF2- 
regulated network should reveal how NRF2 activation 
and suppression can impact processes such as cell homeo- 
stasis, adipogenesis, proliferation, environmental response 
and disease etiology. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Tables 1-6, Supplementary Figures 1-10, 
Supplementary Methods, Supplementary Data set and 
Supplementary References [87-112]. 
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